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Abstract. In this paper we present a new algorithm to compute the Grobner 
basis of an ideal that is invariant under certain permutations of the ring- 
variables. Furthermore, we introduce a second algorithm which is a modifi- 
cation of the modular computation of Grobner bases as introduced by Idrees, 
Pfister, Steidel (cf. HPSllp in the symmetric case. In fact, the algorithm that 
uses the given symmetry, improves the modular calculations in positive char- 
acteristic. In particular, we could, for the first time, compute the Grobner 
basis of the famous ideal of cyclic 9-roots (cf. |BF91| ) over the rationals with 
Singular. Both new algorithms are implemented in Singular. 



Computing the Grobner basis of an ideal is a powerful tool in commutative al- 
gebra, with applications in algebraic geometry and singularity theory. The first 
general algorithm was proposed by B. Buchberger in 1965 (cf. |Bu65| ). Within 
this article we improve the computation of Grobner bases in case that the in- 
put ideal has some special symmetry-character. Consider, for example, the ideal 
/ := (x 2 y 2 — z, xy — 2y + 3z, xy — 2x + 3z) C Q[x, y, z]. Then / does not vary if 
one interchanges the variables x and y, and we say that I is symmetric with respect 
to the permutation x < — > y. In the following we use this property to manipulate 
the ideal by an appropriate linear transformation and apply Buchberger's algorithm 
subsequently. We start in Section [2] by presenting some basic notations, definitions 
and results. In Section [3] we introduce the symmetric Grobner basis algorithm, 
and state a theoretical result that reserves the impact of the symmetry in Propo- 
sition 13.81 Section [4] combines the symmetric algorithm of Section [3] with modular 
methods. Examples and timings are presented in Section [5] 



Let a £ % n :— Sym({l, . . . , n}) be a permutation. The order of a is the minimal 
natural number k G N>o such that o~ k = id, in particular ord(cr) := #((cr)) < oo. 
In order to describe a properly, we make use of the following well-known result 
concerning the representation of permutations. 

Definition 2.1. Let a £ S n be a permutation. Then there exists a natural number 
f9(cr) and a finite disjoint partition {1, . . . , n} = Ufi.^ . . . , e^^} such that 
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with + + = n and < U < n for all i G {1, . . . , i9(cr)}. The cycles 
(ei t i . . . e^/J are up to alignment uniquely determined, and we call this representa- 
tion the cycle decomposition of a. The tuple (h, . . . , is called the cycle type 

Of <T if ll < . . . < 

Note that having the cycle decomposition of a permutation a it holds ord(er) = 
lcm(/i, . . . , From now on we assume that all considered permutations a G S„ 

are given in cycle decomposition. 

Now let K be a field and X := {ari, . . . ,x n } be a set of indctcrminatcs, then 
ffGSn induces a canonical automorphism on the polynomial ring over K in these 
indeterminates, -K^-X"], via tp a : K[X] — > K[X], Xi i — > x a ^y By abuse of notation 
we always write a instead of ip a , i.e. we identify the group S„ as a subgroup of the 
automorphism group Ant(K[X]). 

Definition 2.2. Let / C K\X] be an ideal and a G Aut(K[X]) be an automor- 
phism. Then / is called a-symmetric if a(I) = I. Moreover, let S C Aut(i^[X]) be 
a subgroup then we call / S- symmetric if it is cr-symmetric for all a G S. 

Every subgroup of $„ has only finitely many elements and is therefore finitely 
generated. Hence, let S = (a\, . . . ,a m ) C $„ then an ideal / c K[X] is S- 
symmetric if and only if it is <7i-symmetric for all i G {1, . . . ,m}. In particular, if 
an ideal is cr-symmetric then it is (cr)-symmetric. 

Moreover, given an ideal / C K[X] we can always choose a finite set of poly- 
nomials Fi = {/i,...,/ r } such that / = (Fi). Thus, if / is cr-symmetric with 
cr G Ant(K [X]) we even may assume that cr(F/) = Fi by possibly adding some 
polynomials to F/. 

Example 2.3. The ideal / = (x 2 y 2 - z, xy - 2y + 3z, xy - 2x + 3z) C y, z] is 
obviously cr-symmetric for a — (12)(3) G S3. 

We denote by Mon(X) the set of monomials. Moreover, if > is a monomial 
ordering and / G K[X] a polynomial, then we denote by LC(/) the leading co- 
efficient of /, by LM(/) the leading monomial of /, by LT(/) the leading term 
(leading monomial with leading coefficient) of /, and by tail(/) = / — LT(/) the 
tail of / with respect to the ordering >. In particular, with our notation it holds 
LT(/)=LC(/)-LM(/). 

Convention 2.4. In the following > is a degree ordering, and we always consider 
reduced Grobner bases G, that is ^ G, LM(g) \ LM(/) for any two elements / 7^ g 
in G, and LC(c/) — 1 respectively no monomial of tail(g) is contained in the leading 
ideal of G for any g G G. 

3. Computing Grobner bases using symmetry 

Within this section we describe how to achieve an improvement of the Grobner 
basis computation of a cr-symmetric ideal / G K[X] by using its symmetric prop- 
erty. The basic idea is the construction and usage of an appropriate linear transfor- 
mation t G Aut(K[X)) which "diagonalises" a and respects the cr-symmetry of /. 
It turns out that in many cases the usual Grobner basis computation on the trans- 
formed side is much faster than the computation on the original side. Since the pull 
back of this Grobner basis is in general not a Grobner basis anymore we have to 
add another Grobner basis computation. Nevertheless, this indirection effects an 
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enormous speed up compared to the usual Grobner basis algorithm (cf. Section [5]) . 
We assume that the tupel (a, K) with ord(cr) = k 6 N always satisfies char(X) { k 
and K has a fc-th primitive root of unity 

Remark 3.1. We can always achieve this assumption by possibly adjoining In 
particular, we can swap to K [£/J by working over the field K[a]/<&k(a) where $fc(a) 
is the A;-th cyclotomic polynomial. 

We start by illuminating the basis for the symmetric Grobner basis algorithm 
from a character theoretical point of view and in terms of linear algebra. 

Therefore, we consider the n-dimensional if-subvector space V — (x%, . . . , x n ) K 
of the infinite-dimensional if- vector space ^[-X - ]. Then due to our assumption 
that char(iT) \ #((cr}), character theory guarantees that every representation of 
(a) C S n is a direct sum of irreducible representations (cf. |Se961 Theorem 2]), and 
all irreducible representations of (a) C S„ have degree 1 since (a) C ffi n is an abelian 
group (cf. [Se96, Theorem 9]). In particular, the representation p : (a) — > Aut(V^) 
is diagonalisable, i.e. V = ©" =1 Vi with Vi — (yi) K and p{a)(yi) — £^ ■ yt for some 
Q<Vi<k-l. 

In terms of linear algebra we have the following quite simple proposition which, 
together with its proof, forms the basis of the symmetric Grobner basis algorithm. 

Proposition 3.2. Let a G S„ have cycle type (h, . ■ ■ , £#( CT ))- Then a £ Aut(V) is 
diagonalisable with eigenvalues * m ^' 1 < to < < j < l m — 1}. 

Proof. Let a = (e% . . . e n ) £ $„ with {ei, . . . , e n } — {1, . . . ,n} be an n-cycle. The 
columns of the representation matrix M(cr, X) of a 6 Aut(T^) with respect to 
the iif-basis X = {x\, . . . ,x n } of V are just the permuted unit vectors of K n . 
Hence, M(a,X) is a unitarian matrix and in particular diagonalisable. Moreover, 
let C = tt n — M(a,X) 6 Mat(n x n,K[t]) then the characteristic polynomial of 
a e Aut(V) is 



and {1, £ n , . . . , 1 } are exactly the eigenvalues of a. Now, consider the com- 
binatorial set 



Note that Y is a i-T-basis of V since the coefficients of each y &i are just different 
powers of the primitive root of unity £„, and {x CT j-i( e .) |l<j<n} = Xisa 
if-basis of V itself. 

Let i e {1, n}, then we easily compute that a(y ei ) = -y ei - Consequently, 
yi is the eigenvector corresponding to the eigenvalue , and the representation 
matrix M(a,Y) of a 6 hut(V) with respect to Y is diagonal. 

Now, let a £ S„ have cycle type (1%, . . . , lM a )) with cycle decomposition a = 
°i • " • o^o) and o rn = (e mt i . . . e rri} i m ) for 1 < m < $(cr). Then we have ord(er) = 
k = lcm(/i, . . . , ord(cr m ) = l m , and 




r + sign(a)-(-l) 



n 



<™ + (-l)™- 1 •(-!)"=*"-! 
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is an l m -th. primitive root of unity. We set X m = {x em t , ■ ■ ■ , x em , } such that 
(j m £ Aut(Kn) with V m — (X m ) K , and X = X% U . . . U X$( a ) is a if-basis of 
V = V\ ® . . . © . Analogously to the n-cycle case we obtain the combinatorial 
sets 

lm 



Yn 



IJe 



(i-l)O'-l) 



1< i < L 



of eigenvectors of a m £ Aut(V^ ra ) so that the representation matrix M(a m ,Y m ) of 
er m £ Aut(V m ) with respect to the if-basis Y m of V m is diagonal with eigenvalues 
{l,£j m ,£f > • • • i ^ }• Hence, Y" = Yi U . . . U Y#( CT ) is a if -basis of V, consists of 
eigenvectors of a £ Aut(V), and the representation matrix 

/ M(ai,Yi) \ 
M(o-,y)= •-. eMat(nxri,]f) 

V M(<T^ ia) ,Y^ (a )) J 

of a £ Aut(V) with respect to Y is diagonal with eigenvalues 1 < m < 

#{<r), 0<j<l m ~l} = {£ k/lm)j I 1 < m < < j < l m - 1}. □ 



Remark 3.3. Due to the constructive proof of Proposition 13.21 the eigenvectors of 
a £ Aut(V) can be obtained purely combinatorial which is quite profitable from 
the algorithmic point of view. Hence, let us define the ring homomorphism 



r : K[X] — > K[X] 



ST t('- x )0'-i) 



x e m ,i 1 ► Ue m ,i — / Ah ' U/ <7jT 1 (e m , ( ) 

which maps the ring- variables onto the eigenvectors of a £ Aut(V). Consequently, 
we can define another ring-homomorphism a T induced by the ring- automorphism 
a £ Aut(K[X]) and the linear transformation r £ Aut(K[X]) obtained by the 
commutative diagram 

K[X] — — > K[X] 



■* K [ x ] 



K[X] 

so that a T = rar^ 1 satisfies the property that a T (xi) — £Y* -Xi for suitable exponents 
< Vi < k — 1 and all 1 < i < n. 

Example 3.4. Let a = (12) (3) £ S 3 with ord(er) = 2. Then consider £ 2 = -1 £ Q 
and construct 

T:Q[x,y,z]^Q[x,y,z] 

\0-0 



y 

Z I 



(— 1) • x + (—1) ■ y = x + y, 
(-l) 1 -°-y + (-l) hl -x = y-x, 



as in Remark 13.31 Hence, r is bijective with inverse r 1 defined by 



', T- 1 (Z)=Z. 



GROBNER BASES OF SYMMETRIC IDEALS 



5 



Referring to Remark 1 3. 3 1 ov is induced by o T = tot 1 and thus it holds 



o T {x) = tot (x) = TO 



x— y\ I y — x\ y—x—x+y 



2 J V 2 / 2 

x + y\ ( x + y\ x + y + (y — x) 



°r{y) = tot (y) ( } 

o T (z) = T(TT^ 1 (y) = tct(z) = t(z) = z. 

Notation 3.5. For a better understanding we will index objects that live on the 
transformed side by r. 

As aforementioned respectively proven above, the induced automorphism o T has 
a nice multiplication property on the ring variables what, however, is a priori not a 
sufficient reason for a fast Grobner basis computation. But, in addition, the linear 
transformation r also respects the symmetry of the input ideal. 

Proposition 3.6. If the ideal I C K [X] is o -symmetric, then the transformed ideal 
I T := r(J) G ^[-X] is o T -symmetric. 

Proof. Let / = • • • , /?-)• By definition of o T we obtain for o(fi) — fj that 
o- T «fi)) = r(o-{fi)) = r(/i). Thus, the ideal I T := r(7) - (r(/i), . . . , r(/ r )> is 
cr r -symmetric. □ 

Example 3.7. Let >=>d P be the degree reverse lexicographical ordering, and 
er, r, cr T as in Example 13.41 Now we consider the cr-symmetric ideal 

I = (x 2 y 2 — z, xy — 2y + 3z, xy — 2x + 3z) C Q[x, y, z] 

and obtain that the transformed ideal 

I r := t(I) = (x 4 - 2x 2 y 2 +y A - z, - x 2 + y 2 + 2x - 2y + 3z, 

- x 2 + y 2 - 2x - 2y + 3z) C Q[x, y, z] 

is cr T -symmetric. 

Due to Proposition 13.21 and Proposition ^. 61 we see that transforming the original 
ideal via r still respects some symmetry. In particular, the ideal I T is ov-symmetric 
and applying the automorphism o T on any variable respectively monomial effects 
just a multiplication by a power of a primitive root of unity. The advantage of 
this circumstance is the fact that the symmetry propagates during the process of 
computing a Grobner basis of I T which influences the performance in a positive 
way. More precisely, the following proposition holds. 

Proposition 3.8. Let I T be o T -symmetric, then a Grobner basis G T of I T satisfies 
a T (g T ) — £, 9 ' r • g T for all g T G G T with suitable < v 9t < k — 1. 

Proof. Let I T = (Fj t ) and /, g G Fj t . Due to the property of o T we can define 
X a := LM(/) = LM(o- T (/)), X? := lM{g) = LM(o T (g)), X~< := lcm(X a ,X^) and 
o T {X a ) = ^X a , o T {X?) = QX?, o T (X~<) = ftXi for suitable u ai v p ,u^ G 
{0,...,fc-l}. Then 

spoly(/, 5 ) = ■ f - ^gj • X^ ■ g 



^Degree reverse lexicographical ordering: Let X a ,X^ £ Mon(A'). X a >dp X@ :<^=> 
deg(X«) > deg(XP) or (deg(X a ) = deg( X^) an d 31 < i < n : a n = ,<*i_i = ft-i.a* < 

ft), where deg(X a ) = «i + . . . + £*„; cf. IGP07| . 
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and again by the property of a T it holds LC(ov(/)) = Cfc" ' LC(/) respectively 
LC(a T (g))=C -^C(g). Thus 

ov(spoly(/,ff)) = a r pr-«) • <r T (/) - • a T (X^) • a r (.g) 

= £,-„„ x ^ a aAf) _ LCO) x7 _, 

t v~,-v a ( vl - a , j.s LC(<7 T (/)) / n 

= C ^ • s POly(CTr(/),CT r (g)). 

Moreover, there are a^,r G such that spoly(/,<?) = X^gf, + r - Due to 

the above computation it follows 

spoly(CT T (/),o- T (5)) = C ' °V X! a >^ + r = X! b h h + £k a ^tW, 

for suitable 6^ = ^ T • a^-i^ G ^[^] since I T respectively -Fj T is cr T -symnietric 
and consequently 

NF ( spoly(o> (/) ,a T {g)),F lT ) = ' ^ (NF ( spoly(/, g), i*j T )) . 

This property implies that the reduced Grobner basis G T = . . . , g T s } of I T 
satisfies a T (gJ) = £^ lJ ■ gj for suitable i, j G {1, . . . , s} and Vij G {0, . . . , k — 1}. 
Moreover, it follows LM(gJ) = LM(cr T (g l r )) = LM(<?J), but since G r is reduced we 
conclude gj — gj. Hence, we have a T (gf) = ■ gj for all i G {1, . . . , s} with 
suitable ^ G {0, . . . , k — 1}. □ 

Example 3.9. Let I T = (x 4 — 2x 2 y 2 + y 4 — 16z, x 2 — y 2 — 4a; + Ay + 12z, x 2 — y 2 — 
Ax — Ay + 12z) C Q[a;, y, z] and cr T G Aut(Q[a;, y, z]) as obtained in Example 13.71 
Then I T is ov-symmetric, and its Grobner basis 

G T = {x, 12yz - 9z 2 - 8y + 13z, y 2 - 2y + 3z, 81z 3 + 36z 2 - 56y + 115 z} 

satisfies cr T (g T ) = (— l) 1 ^ • g T for suitable v 9t G {1,2} and all g T G G T . Now, the 
reverse transformation of G T yields the set 

r _1 (G r ) = {\ x — \y, 6a;z + 6yz — 9z 2 — Ax — Ay + 13z, 

\x 2 + \xy + \y 2 - x - y + 3z, 81z 3 + 36z 2 - 28a; - 28y + lWz}. 

Obviously, just pulling back the Grobner basis G T via r _1 does not lead to a 
Grobner basis of the input ideal /. Thus, we have to compute a Grobner basis of 
the ideal (t _1 (G t )) as well. Nevertheless, the advantage of this computation is 
the fact that the achieved property as described in Proposition 13.81 is respected by 
applying r^ 1 on G T . More precisely, the following proposition holds. 

Proposition 3.10. a(g) = ^ a ■ g for all g G r _1 (G T ) and suitable < v g < k — 1. 

Proof. Let g G r _1 (G r ), i.e. there is an g T G G T such that g = r^ 1 (g T ). Then due 
to Proposition 13.21 and Proposition 13.81 we have 

r(a(g)) = r{a{T-\g T ))) = a T (g T ) = • g T 
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for some v 9r £ {0, . . . , k — 1}. Hence, we obtain 

a(g) = T-\r{a{g))) = r" 1 • 3.) = C T ' ^(flr) = ^ ' 5- 
This proves the proposition. □ 

Example 3.11. Let a = (12) (3) £ S 3 with ord(cr) = 2, £ 2 = -1 £ Q and T- l (G T ) 
as obtained in Example 13.91 Then we compute 

a{\x - \y) = -{\x - \y), 

a(6xz + 6yz — 9z 2 — Ax — 4y + 13 z) = 6xz + 6yz — 9z 2 — Ax — Ay + 13z. 

a (i x2 + \ x v + \v 2 - x - v + 3z ) = \ x2 + h x V + W ~ x ~ V + 3z ' 

cr(81z 3 + 36z 2 - 28x - 28y + 115 z) = 81z 3 + 36z 2 - 28a; - 28y + 115z, 

as claimed in Proposition 13. 101 

The following diagram summarizes and illustrates our way of improving the 
computation of a Grobner basis G of a c-symmetric ideal /. 

(I, a) — — > (l T ,a T ) 

std 

(t-HGt)) G T 
| std 

G 

Note that the linear transformation r is defined in Remark 13.31 and the proce- 
dure std is implemented in Singular and computes a Grobner basis respectively 
standard basis of the input. 

Algorithm [T] computes the Grobner basis of a er-symmetric ideal 70 

Algorithm 1 Symmetric Grobner Basis Computation (symmStd) 

Input: I C if[A], a £ S„ such that I is er-symmetric. 
Output: G C K[X] the Grobner basis of I. 

1: k := ord(cr); 

2: if (k mod char(i^) = 0) then 

3: print Warning, algorithm is not applicable. 

4: return 0; 

5: if (k = 2 or (char(AT) - 1) mod k = 0) then 
6: compute £ K; 
7: else 

8: iT = K[a]/$ k (a); 
9: £ fe := a: 

10: compute r £ Aut(K[X]); 
11: compute a Grobner basis G T of 7 T = r(J); 
12: compute a Grobner basis G of (t _1 (G t )); 
13: return G; 



'The corresponding procedures are implemented in Singular in the library symodstd. lib. 
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Theorem 3.12. Algorithm^ terminates and is correct, i.e. the output G is a 
Grobner basis of the input I . 

Proof. Termination is clear and for proving correctness it suffices to show that 
I = (t~ x (G t )) since G is by definition a Grobner basis of (r _1 (G r )). Let / 6 
/ and G T = {g[ , . . . , gj}. Then r(/) 6 t{I) = I T and consequently there are 
oi, . . . , a s € K [X] such that r(/) = a * ' ffl since G T is a Grobner basis of I T . 

Hence, we obtain 

/ = r-\r(f)) = J2 T-H<H) ■ r-\9j) e (r-\G T )) . 
»=i 

For the other inclusion let g G (r _1 (G T )). It follows that r(g) G (G T ) = I T = t(I) 
and moreover g (z I since r is an automorphism. □ 

For illustration of Algorithm [T] we combine all previous examples. 

Example 3.13. Again, let / = (x 2 y 2 — z, xy — 2y + 3z, xy — 2x + 3z) C y, z] 
and a = (12)(3) £ S 3 . Referring to Examples El E2I EU and [3TT] we already 
obtained 

I T := r(J) = (a; 4 - 2x 2 y 2 + y 4 - z, - x 2 + y 2 + 2x - 2y + 3z, 

-x 2 + y 2 -2x-2y + 3z), 

and its Grobner basis 

G T = {x, \2yz - 9z 2 - 8y + I3z, y 2 - 2y + 3z, 81z 3 + 36z 2 - 56y + U5z} 
with 

r _1 (G T ) = {\x - \y, 6xz + 6yz - 9z 2 - Ax - Ay + 13z, 

ia; 2 + \xy + \y 2 - x - y + 3z, 81z 3 + 36z 2 - 28x - 28y + 115z}. 

Finally, we compute 

G = {x-y, \2yz ~ 9z 2 - 8y + 13z, y 2 - 2y + 3z, 81z 3 + 36z 2 - 56y + 115z}, 

the Grobner basis of (t _1 (G t )) respectively I. 

Remark 3.14. It is remarkable that we even achieve an advancement although we 
have to compute Grobner bases internally twice on modified input ideals via std 
(cf. Section EJ. 



4. Computing Grobner bases using symmetry and modular methods 

When applying Algorithm [T] on <7-symmetric ideals over the rationals such that 
ord(cr) = k > 2 we need to swap to Q[a]/$fe(a) as explained in Remark 13.11 How- 
ever, over fields K of positive characteristic such that k | (char(A') — 1) this can be 
omitted. Consequently, we use modular methods to improve Algorithm [T] applied 
on er-symmetric ideals in the polynomial ring over the rationals. More precisely, we 
improve the modular Grobner basis algorithm as introduced by Arnold (cf. [A03 ) 
and Idrees, Pfister, Steidel (cf. [IPSllj l. 
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Algorithm [2] combines the algorithms symmStd (cf. Algorithm!!]) and modStd (cf. 
[IPS11I Algorithm 1])@ 



Algorithm 2 Symmetric Modular Grobner Basis Computation (syModStd) 

Input: / C Q[X] and a G S n such that I is cr-symmetric. 
Output: G C Q[X] the Grobner basis of I. 

1: k := ord(cr); 

2: choose P, a list of random primes such that k \ (p — 1) for all peP; 
3: GP 0; 

4: loop 

5: for p E P do 

6: G p = symmStd(/p, cr); 

7: GP = GPU{G P }; 

8: (GP P) = deleteUnluckyPrimesSB(GP, P); 

9: lift (GP, P) to G C Q[X] by applying Chinese remainder and Farey rational 
map; 

10: if pTestSB(I, G, P) then 

li: if J C (G) then 

12: if G is a standard basis of (G) then 

13: return G; 

14: enlarge P; 



Remark 4.1. The essential differences compared to the algorithm modStd are the 
following: 

(1) The choice of the prime list P has to be restricted. Every considered prime 
number p G P has to satisfy the condition k \ (jp — 1) in order to assure that 
the coefficient field ¥ p has a fc th primitive root of unity. 

(2) The modular Grobner bases G p are computed via symmStd instead of std. 

Similar to modStd we can parallelize Algorithm[5]by computing the modular Grobner 
bases G p respectively performing the final tests in parallel. 

Theorem 4.2. Algorithm [H terminates and is correct, i.e. the output G is a 
Grobner basis of the input I . 

Proof. Termination is clear and correctness follows directly from Theorem 13.121 and 
|IPsill Theorem 2.4]. □ 

5. Examples and timings 

In this section we provide examples on which we time the new algorithms 
symmStd (cf. Section [3) respectively syModStd (cf. Section 2]) and its paralleliza- 
tion as opposed to the former algorithms std respectively modStd implemented in 
Singular. Timings are conducted by using Singular 3-1-3 on an AMD Opteron 
6174 machine with 48 CPUs, 2.2 GHz, and 128 GB of RAM running the Gentoo 
Linux operating system. 



^The corresponding procedures are implemented in Singular in the library symodstd.lib. 
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Example 5.1 (Cyclic n-roots (cf. |Bj85| , |Bj90| , [BF91] )). The task to compute a 
Grobner basis of the ideal in Q[X] = Q[xi, . . . ,x n ] corresponding to the following 
system of polynomial equations 

xx + . . . + x n = 
X1X2 + X2X3 + . . . + X n -iX n + x n xx = 

XiX 2 ■ • ■ X n -l + X 2 X 3 ■ ■ ■ X n + . . . + X n ^iX n ■ ■ • X„_ 3 + X n Xx • • • X n - 2 = 

xx ■ ■ ■ x n — 1 = 

has become a benchmark problem for Grobner basis techniques. We call this ideal 
cyclic(n), and its variety cyclic n-roots (cf. |Bj85]). The origin of the problem is 
related to Fourier analysis (cf. |Bj85| , |Bj90 ). Obviously, the ideal cyclic(n) is by 
definition symmetric with respect to the n-cycle a n = (1 . . . n) such that we can 
apply the algorithms symmStd and syModStd. 

Until the end of 2009, Singular was able to compute a Grobner basis of cyclic(n) 
for n < 8. In April 2010, we could for the first time compute a Grobner basis of 
cyclic(9) via a prototype of syModStd by using the 32-bit version of Singular 3-1-1 
on an Intel® Xeon® X5460 machine with 4 CPUs, 3.16 GHz each, and 64 GB of 
RAM under the Gentoo Linux operating system within 23 days. 

Table[T]summarizes the present timings for computing a Grobner basis of cyclic(n) 
for n — 7, 8,9 with different numbers of cores t denoted by modStd(£) respectively 
syModStd(^), and different permutations a where again k — ord(er) denotes the 
order of a. 



n 


k 


symmStd 


modStd(l) 


modStd(30) 


syModStd(l) 


syModStd(30) 


7 


7 


409.394 


106 


38 


165 


50 


8 


2 




6.816 


973 


5.196 


811 


8 


4 




6.816 


973 


3.120 


620 


8 


8 




6.816 


973 


4.454 


788 


9 


9 




9.935.103 


475.981 


2.790.303 


207.681 



Table 1 . Total running times in seconds for computing a Grobner 
basis of cyclic(n) for n = 7,8,9 via symmStd, modStd(^), and 
syModStd(^) for £ = 1,30, using different permutations of order 
k. The symbol "-" indicates out of memory failures. 



In these examples we used the permutation (1234567) G S7 for n = 7, the 
permutations (18)(27)(36)(45), (1753)(2864), (12345678) 6 $s for n = 8, and the 
permutation (147) (258) (369) G $9 for n — 9. Note that the timings obtained by 
the modular versions are dependent on the used permutation and especially on its 
order. In particular, a higher order fc, that is a higher symmetry, speeds up the 
Grobner basis computation on the transformed side but in contrast slows down the 
application of the linear transformation r (cf. Remark I3.3[) since the support of a 
ring- variable's image depends on the order k of the permutation. This circumstance 
justifies that applying the symmetric modular algorithm for computing a Grobner 
basis of cyclic(8) is most performant when using the permutation (1753)(2864) G S>s 
of order 4. 
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Example 5.2 (100 Swiss Francs Problem (cf. jZJGllj . [St08] )). B. Sturmfels offered 
a cash prize of 100 Swiss Francs for the resolution of a very specific conjecture in 
the Nachdiplomsvorlesung (postgraduate course) which he held at ETH Zurich in 
the summer of 2005. Based on a concrete biological example proposed in [PS051 
Example 1.16] the problem arised to maximize the likelihood function 

-40 

hp)=[Upu) [Up*] ■ ( E*« 

over all (positive) 4x4- matrices P = (py)i<*,i<4 °f rank at most two. Due to 
numerical experiments by applying an expectation-maximization algorithm (EM 
algorithm), B. Sturmfels conjectured that the matrix 



40 



/ 3 3 2 2 \ 

3 3 2 2 

2 2 3 3 

\ 2 2 3 3 ) 



is a global maximum of the likelihood function L(P) (cf. [St08] ). 

The conjecture is positively confirmed in [ZJGllj . In their approach via Grobner 
bases (cf. [ZJG11, Section 2.3]) it is necessary to compute the Grobner basis of the 
ideal J defined by 



4 4 



h,^2ai,^2bi, f -l,..., f 4 , g-i, gij C Q[ ai , ...a 4 ,b 1: . 



=1 i=l 



with 

/i = E (bj-il + aibj-llil + aibk)] + h ■ JJ (1 + ath) 
j=i \ k^j J fc=i 

4 / \ 

9i = / A a j ■ (1 + a ibi) ■ JJ(1 + a k bi) + flj • JJ(1 + a k bi) 
j=i \ k^j J k=l 

for 1 < i < 4, and 

J = I + (1 -ttoi) C Q[a 1 ,...,a i ,b 1 ,...,b 4 ,u] 

with respect to an elimination ordering on the variable u. In a first approach 
we therefore applied modStd using the lexicographical ordering >i p respectively 
the block ordering (>dp(8)i to eliminate the variable u. It turned out that 

both variants are comparably slow so that we used in a second approach the de- 
gree reverse lexicographical ordering >d P , and applied the FGLM-algorithm (cf. 
FGLM93 ) subsequently to obtain a Grobner basis with respect to the block or- 
dering (>dp(8), >; P (i))- Since the ideal J C Q[a±, . . . , 04, b\, . . . , 64, u] is symmetric 
with respect to the permutation (34) (78) € S9 we could moreover apply the algo- 
rithm syModStd. The timings for the computations in Singular are summarized 
in Tabled 
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Method 


Running Time 


modStd[>! p ] 


39919 


modStd[(> dp (8), >i p (i))] 


515 


syModStd[(> dp(8) ,> ;p(1) )] 


356 


modStd[> dp ] - fglm[(> dp ( 8 ), > ;p{ i))] 


375 


syModStd[> dp ] - fglm[(> dp(8) , > ip(1) )] 


284 



Table 2. Total running times in seconds for computing the 
Grobner basis of J C Q[ai, . . . , 04, b\, . . . , &4, u] with respect to 
an elimination ordering on the variable u via different methods. 



Example 5.3 (Inverse Galois Problem (cf. |Mal94j , |Mat87j . [MM99] )). A major 
topic in algebraic number theory is the inverse Galois problem over a field K, i.e. 
the question whether any finite group G is the Galois group of a Galois extension 
of K. The most interesting case is K = Q which is still open in general. In 
contrast, the problem is known to be true for K being a rational function field in 
one variable t over an algebraically closed field of characteristic zero. In particular, 
it is true for K = C(t), and in this case it is solved via geometric field extensions 
(see for example |MM99( I, §1]). Moreover, the same strategy applies to finite 
field extensions of Q(t) ramified only over {0,1, 00} (see for example [MM99, I, 
§5]). In this situation, for any triple a — (<7i, a<i, 03) of elements generating a 
transitive subgroup G = (a%, 02, a 3) C $„ with <j\<JiOz = 1 there exists a certain 
field extension K a /Q(t) of degree n, unramified outside {0, 1, 00}, and whose Galois 
group is isomorphic to G. In fact, any such extension K a /Q(t) is already defined 
over a number field k a = Q(a a ), the so-called field of definition of K a /Q(t), so that 
there exists a further field extension K a lk a (£) which also has G as its Galois group. 
The degree [k a : Q] is bounded from above by group theoretical information (see 
for example |Mal941 Proposition A]). In order to construct the extension K a /Q(t) 
it is necessary to solve a system of polynomial equations (see MM99, I, §9]). In 
case that, for example, o\ and 02 have the same cycle type, the defining ideal 
is symmetric with respect to a permutation of order 2 so that we can apply the 
algorithms symmStd and syModStd to compute a Grobner basis of this system. 
In addition, choosing an elimination ordering for the last ring variable, the last 
polynomial / of the Grobner basis of this system of non-linear equations generates 
the field of definition k a = Q(a a ) insofar that a a is a zero of /. The irreducible 
factors of / together with group theoretical information yield restrictions on [k a : 
Q]. In case that [k a : Q] = 1, the given group G can even be realized over Q(t), 
and therefore also over Q by Hilbert's irreducibility theorem. 

In 1994, G. Malle collected computational data on several k a of degree [k a : Q] < 
13 (cf. |Mal94j ) with the intention to observe regularities and hints to decrease the 
group theoretical bound. Table [3] lists further examples in the spirit of this article 
and which could not be computed at that time. 

Note that all ideals belonging to the examples listed in Table[3]are zero-dimensio- 
nal such that we can compute a Grobner basis with respect to the degree reverse 
lexicographical ordering, and obtain a lexicographical Grobner basis by applying 
the FGLM-algorithm (cf. |FGLM93j ) subsequently. 
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n 


G 


c a 


[** : Q] 


symmStd 


modStd 


syModStd 


7 


A 7 


4.2 - 4.2 - 4.2 


12 


24 


19 


19 


9 




4.2^ - 4.2^ — 5.3 


34 


5 


15 


13 


10 


Aio 


5.2^ - 5.2^ - 7 


37 


977 


186 


107 


11 


An 


4.2-9-9 


12 


17 


2 


1 


11 


An 


4.2-4^.3-4^.3 


8 


34 


3 


2 


11 


An 


4.2 - 5.3^ - 5.3^ 


8 


15 


3 


2 


11 


An 


5 - 8.2 - 8.2 


11 


265 


17 


8 


11 


An 


5 - 6.4 - 6.4 


11 


339 


20 


10 


11 


An 


5 - 7.3 - 7.3 


11 


292 


16 


8 


11 


Sn 


7-4.3^ -4.3* 


26 


631245 


2506 


1493 


11 


Sn 


7-4.3.2^ -4.3.2 2 


29 




3039 


1979 


11 


Sn 


7 - 7.2 - 7.2 


29 




1414 


702 


11 


ffill 


7-5.4-5.4 


26 




1738 


899 



Table 3. Total running times in seconds for computing the defin- 
ing Grobner basis of the field extension K a /Q(t) having group G 
and conjugacy class triple C a (cf. |Mal94] ) via symmStd, modStd, 
and syModStd. Here, C a is a class of G containing elements of the 
given cycle type. The symbol "-" indicates out of memory failures. 



Conclusion 5.4. In all considered examples the symmetric modular version of the 
Grobner basis algorithm is the most performant one, and is, hence, a quite powerful 
tool if the input ideal is symmetric with respect to some permutation of the ring- 
variables. 
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